Search Results

The current trend of downsizing used in gasoline engines, while reducing fuel consumption and CO2 emissions, imposes severe thermal loads inside the combustion chamber. These critical thermodynamic conditions lead to the possible auto-ignition (AI) of fresh gases hot-spots around Top-Dead-Center (TDC). At this very moment where the surface to volume ratio is high, wall heat transfer influences the temperature field inside the combustion chamber. The use of a realistic wall temperature distribution becomes important in the case of a downsized engine where fresh gases hot spots found near high temperature walls can initiate auto-ignition. This paper presents a comprehensive numerical methodology for an accurately prediction of thermodynamic conditions inside the combustion chamber based on Conjugate Heat Transfer (CHT).

Actual combustion strategies in internal combustion engines rely on fast and accurate injection systems to be successful. One of the injector designs that has shown good performance over the past years is the direct-acting piezoelectric. This system allows precise control of the injector needle position and hence the injected mass flow rate. Therefore, understanding how nozzle flow characteristics change as function of needle dynamics helps to choose the best lift law in terms of delivered fuel for a determined combustion strategy. Computational fluid dynamics is a useful tool for this task. In this work, nozzle flow of a prototype direct-acting piezoelectric has been simulated by using CONVERGE. Unsteady Reynolds-Averaged Navier-Stokes approach is used to take into account the turbulence. Results are compared with experiments in terms of mass flow rate. The nozzle geometry and needle lift profiles were obtained by means of X-rays in previous works.

The 3-Zones Extended Coherent Flame Model (ECFM3Z) and the Tabulated Kinetics for Ignition (TKI) auto-ignition model are widely used for RANS simulations of reactive flows in Diesel engines. ECFM3Z accounts for the turbulent mixing between one zone that contains compressed air and EGR and another zone that contains evaporated fuel. These zones mix to form a reactive zone where combustion occurs. In this mixing zone TKI is applied to predict the auto-ignition event, including the ignition delay time and the heat release rate. Because it is tabulated, TKI can model complex fuels over a wide range of engine thermodynamic conditions. However, the ECFM3Z/TKI combustion modeling approach requires an efficient predictive spray injection calculation. In a Diesel direct injection engine, the turbulent mixing and spray atomization are mainly driven by the liquid/gas coupling phenomenon that occurs at moving liquid/gas interfaces.

In internal combustion engine computational fluid dynamics (CFD) simulations, uncertainties arise from various sources, such as estimates of model parameters, experimental boundary conditions, estimates of chemical kinetic rates, etc. These uncertainties propagate through the model and may result in discrepancies compared to experimental measurements. The relative importance of the various sources of uncertainty can be quantified by performing a sensitivity analysis. In this work, global sensitivity analysis (GSA) was applied to engine CFD simulations of a low-temperature combustion concept called gasoline compression ignition, to understand the influence of experimental measurement uncertainties from various sources on specific targets of interest-spray penetration, ignition timing, combustion phasing, combustion duration, and emissions. The sensitivity of these targets was evaluated with respect to imposed uncertainties in experimental boundary conditions and fuel properties.

Reynolds-averaged Navier-Stokes (RANS) modeling is expected to deliver an ensemble-averaged result for the majority of turbulent flows. This could lead to the conclusion that multi-cycle internal combustion engine (ICE) simulations performed using RANS must exhibit a converging numerical solution after a certain number of consecutive cycles. However, for some engine configurations unsteady RANS simulations are not guaranteed to deliver an ensemble-averaged result. In this paper it is shown that, when using RANS modeling to simulate multiple engine cycles, the cycle-to-cycle variations (CCV) generated from different initial conditions at each cycle are not damped out even after a large number of cycles. A single-cylinder GDI research engine is simulated using RANS modeling and the numerical results for 20 consecutive engine cycles are evaluated for two specific operating conditions.

The primary strength of large eddy simulation (LES) is in directly resolving the instantaneous large-scale flow features which can then be used to study critical flame properties such as ignition, extinction, flame propagation and lift-off. However, validation of the LES results with experimental or direct numerical simulation (DNS) datasets requires the determination of statistically-averaged quantities. This is typically done by performing multiple realizations of LES and performing a statistical averaging among this sample. In this study, LES of n-dodecane spray flame is performed using a well-mixed turbulent combustion model along with a dynamic structure subgrid model. A high-resolution mesh is employed with a cell size of 62.5 microns in the entire spray and combustion regions. The computational cost of each calculation was in the order of 3 weeks on 200 processors with a peak cell count of about 22 million at 1 ms.

The primary objective of the research discussed here was to compare the commercial computational fluid dynamics (CFD) software, CONVERGE, and a prevalent open-source code, OpenFOAM, with regard to their ability to predict spray and combustion characteristics. The high-fidelity data were obtained from the Engine Combustion Network (ECN) at Sandia National Laboratory in a constant-volume combustion vessel under well-defined, controlled conditions. The experiments and simulations were performed by using two diesel surrogate fuels (i.e., n-heptane and n-dodecane) under both non-reacting and reacting conditions. Specifically, ECN data on spray penetration, liquid length, vapor penetration, mixture fraction, ignition delay, and flame lift-off length (LOL) were used to validate both codes. Results indicate that both codes can predict the above experimental characteristics very well.

It is well known that in-nozzle flow behavior can significantly influence the near-nozzle spray formation and mixing that in turn affect engine performance and emissions. This in-nozzle flow behavior can, in turn, be significantly influenced by fuel properties. The goal of this study is to characterize the behavior of two different fuels, namely, a straight-run naphtha that has an anti-knock index of 58 (denoted as “Full-Range Naphtha”) and n-dodecane, in a simulated multi-hole common-rail diesel fuel injector. Simulations were carried out using a fully compressible multi-phase flow representation based on the mixture model assumption with the Volume of Fluid method. Our previous studies have shown that the characteristics of internal and near-nozzle flow are strongly related to needle motion in both the along- and off-axis directions.

Large-eddy Simulations (LES) have been carried out to investigate spray variability and its effect on cycle-to-cycle flow variability in a direct-injection spark-ignition (DISI) engine under non-reacting conditions. Initial simulations were performed of an injector in a constant volume spray chamber to validate the simulation spray set-up. Comparisons showed good agreement in global spray measures such as the penetration. Local mixing data and shot-to-shot variability were also compared using Rayleigh-scattering images and probability contours. The simulations were found to reasonably match the local mixing data and shot-to-shot variability using a random-seed perturbation methodology. After validation, the same spray set-up with only minor changes was used to simulate the same injector in an optically accessible DISI engine. Particle Image Velocimetry (PIV) measurements were used to quantify the flow velocity in a horizontal plane intersecting the spark plug gap.

Traditional Lagrangian spray modeling approaches for internal combustion engines are highly grid-dependent due to insufficient resolution in the near nozzle region. This is primarily because of inherent restrictions of volume fraction with the Lagrangian assumption together with high computational costs associated with small grid sizes. A state-of-the-art grid-convergent spray modeling approach was recently developed and implemented by Senecal et al., (ASME-ICEF2012-92043) in the CONVERGE software. The key features of the methodology include Adaptive Mesh Refinement (AMR), advanced liquid-gas momentum coupling, and improved distribution of the liquid phase, which enables use of cell sizes smaller than the nozzle diameter. This modeling approach was rigorously validated against non-evaporating, evaporating, and reacting data from the literature.

Using gasoline and diesel simultaneously in a dual-fuel combustion system has shown effective benefits in terms of both brake thermal efficiency and exhaust emissions. In this study, the dual-fuel approach is applied to a light-duty spark ignition (SI) gasoline direct injection (GDI) engine. Three combustion modes are proposed based on the engine load, diesel micro-pilot (DMP) combustion at high load, SI combustion at low load, and diesel assisted spark-ignition (DASI) combustion in the transition zone. Major focus is put on the DMP mode, where the diesel fuel acts as an enhancer for ignition and combustion of the mixture of gasoline, air, and recirculated exhaust gas. Computational fluid dynamics (CFD) is used to simulate the dual-fuel combustion with the final goal of supporting the comprehensive optimization of the main engine parameters.

Numerical modeling of fuel injection in internal combustion engines in a Lagrangian framework requires the use of a spray-wall interaction sub-model to correctly assess the effects associated with spray impingement. The spray impingement dynamics may influence the air-fuel mixing and result in increased hydrocarbon and particulate matter emissions. One component of a spray-wall interaction model is the splashed mass fraction, i.e. the amount of mass that is ejected upon impingement. Many existing models are based on relatively large droplets (mm size), while diesel and gasoline sprays are expected to be of micron size before splashing under high pressure conditions. It is challenging to experimentally distinguish pre- from post-impinged spray droplets, leading to difficulty in model validation.

Combustion cycle-to-cycle variation (CCV) of Spark-Ignition (SI) engines can be influenced by the cyclic variations in charge motion, trapped mass and mixture composition inside the cylinder. A high CCV leads to misfire or knock, limiting the engine’s operating regime. To understand the mechanism of the effect of flow field and mixture compositions on CCV, the present numerical work was performed in a single cylinder Direct Injection Spark-Ignition (DISI) engine. A large eddy simulation (LES) approach coupled with the G-equation combustion model was developed to capture the CCV by accurately resolving the turbulent flow field spatially and temporally. Further, the ignition process was modeled by sourcing energy during the breakdown and arc phases with a line-shape ignition model which could move with the local flow. Detailed chemistry was solved both inside and outside the flame front. A compact 48-species 152-reactions primary reference fuel (PRF) reduced mechanism was used.

Low-temperature heat release (LTHR) in spark-ignited internal combustion engines is a critical step toward the occurrence of auto-ignition, which can lead to an undesirable phenomenon known as engine knock. Hence, correct predictions of LTHR are of utmost importance to improve the understanding of knock and enable techniques aimed at controlling it. While LTHR is typically obscured by the deflagration following the spark ignition, extremely late ignition timings can lead to LTHR occurrence prior to the spark, i.e., pre-spark heat release (PSHR). In this research, PSHR in a boosted direct-injection SI engine was numerically investigated using three-dimensional computational fluid dynamics (CFD). A hybrid approach was used, based on the G-equation model for representing the turbulent flame front and the multi-zone well-stirred reactor model for tracking the chemical reactions within the unburnt region.

In this last decade, non-destructive X-ray measurement techniques have provided unique insights into the internal surface and flow characteristics of automotive injectors. This has in turn contributed to enhancing the accuracy of Computational Fluid Dynamics (CFD) models of these critical injection system components. By employing realistic injector geometries in CFD simulations, designers and modelers have identified ways to modify the injectors’ design to improve their performance. In recent work, the authors investigated the occurrence of cavitation in a heavy-duty multi-hole diesel injector operating with a high-volatility gasoline-like fuel for gasoline compression ignition applications. They proposed a comprehensive numerical study in which the original diesel injector design would be modified with the goal of suppressing the in-nozzle cavitation that occurs when gasoline fuels are used.

A wide variety of spray models and their associated sub-models exist to assist with numerical spray development studies in the many applicable areas viz., turbines, internal combustion engines etc. The accuracy of a simulation when compared to the experiments varies, as these models chosen are varied. Also, the computational grid plays a crucial role in model correctness; a grid-converged CFD study is more valuable and assists in proper validation at later stages. Of primary relevance to this paper are the combustion models for a grid-converged Lagrangian spray modeling scenario. CONVERGE CFD code is used for simulation of split injection diesel (n-heptane) sprays and a structured methodology, using RNG k-ε turbulence model, is followed to obtain a grid-converged solution for the key Computational Fluid Dynamics (CFD) parameters viz., grid size, injected parcels and spray break-up time constant.

For several years there has been a great deal of effort made in researching ways to run a compression ignition engine with simultaneously high efficiency and low emissions. Recently much of this focus has been dedicated to using gasoline-like fuels that are more volatile and less reactive than conventional diesel fuel to allow the combustion to be more premixed. One of the key challenges to using fuels with such properties in a compression ignition engine is stable engine operation at low loads. This paper provides an analysis of how stable gasoline compression ignition (GCI) engine operation was achieved down to idle speed and load on a multi-cylinder compression ignition engine using only 87 anti-knock index (AKI) gasoline. The variables explored to extend stable engine operation to idle included: uncooled exhaust gas recirculation (EGR), injection timing, injection pressure, and injector nozzle geometry.

Global Sensitivity Analysis (GSA) is conducted for a diesel engine simulation to understand the sensitivities of various modeling constants and boundary conditions in a global manner with regards to multi-target functions such as liquid length, ignition delays, combustion phasing, and emissions. The traditional local sensitivity analysis approach, which involves sequential perturbation of model constants, does not provide a complete picture since all the parameters can be uncertain. However, this approach has been studied extensively and is advantageous from a computational point of view. The GSA simultaneously incorporates the uncertainty information for all the relevant boundary conditions, modeling constants, and other simulation parameters. A global analysis is particularly useful to address the important parameters in a model where the response of the targets to the values of the variables is highly non-linear.

Combustion is governed by only two phenomena: chemical reactions and mixing (i.e., transport of energy, species, and momentum). A Reynolds Averaged Navier-Stokes (RANS) turbulence model is commonly employed to account for the enhanced mixing due to the presence of turbulence in fluid flow. A RANS turbulence model enhances mixing by introducing a turbulent viscosity. The addition of a turbulent viscosity not only enhances mixing but it also eliminates smaller scales in the CFD simulation. Even though the turbulent viscosity eliminates smaller scales, it is common for RANS engine combustion simulations to be under-resolved. The lack of sufficient mesh resolution to resolve the remaining scales in a RANS combustion simulation may result in a significant sub-grid term that needs to be modeled. In the context of combustion simulation, it is shown that frequently this sub-grid term is significantly more important than Turbulent Chemistry Interaction terms (TCI).

The occurrence of cavitation inside injectors is generally undesirable since it can cause material erosion and result in deviations from the expected injector performance. Previous numerical work employing an injector geometry measured with x-ray diagnostics and operating with a high-volatility straight-run gasoline (SRG) has shown that: (1) most of the cavitation is generally observed at low needle lifts, (2) needle motion is responsible for asymmetric structures in the internal flow as well as large pressure and velocity gradients that trigger phase transition at the orifice inlets, and (3) cavitation affects the injector discharge coefficient and distribution of injected fuel. To explore the potential for material damage within the injector orifices due to cavitation cloud collapse, the cavitation-induced erosion risk assessment (CIERA) tool has been applied for the first time to the realistic geometry of a heavy-duty injector using the CONVERGE software.